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■ S ■ 1 Introduction 

■ 

I— I' The main purpose of this paper is to discuss a common property of certain classes of quasilinear 

I . differential equations, namely: a Lie point symmetry of the considered equation is a Noether 

I symmetry if and only if the equation parameters assume critical values. 

■ By 'Noether symmetry' we mean a variational or a divergence symmetry. Further, as it is well 
(<~^ . known, the so-called critical exponent is found as the critical power for embedding theorems. It 
T-H I is also related to some numbers dividing the existence and nonexistence cases for the solutions of 

' differential equations, in particular, of semilinear differential equations with power nonlinearities 

Q ■ involving the Laplace operator. The above stated property traces a connection between these 

i two notions: the Noether symmetries and the 'criticality' of the equation. 

I We shall review the results obtained in [H El El El ? where that property was established for 

^ ■ various differential equations. To the author's knowledge, the 'enigmatic' relation between the 

^ , Noether symmetries and the critical exponents has not previously been emphasized. 

Our research was inspired by the paper of Clement, de Figueiredo and Mitidieri [Jj in which 

■ these authors introduced the following class of quasilinear ordinary differential equations: 



a 1/1/3 



IX \y 



y')' = xy{y), (1) 



where a, /3 and 7 are real numbers, x > 0, y = y{x) and / is a non-negative function. The 
relations between these parameters will be specified below. We point out that for specific values 
of a, /9, 7 it contains well-known important differential equations briefly listed in Section 2. 
In particular, it contains the radial forms of partial differential equations involving Laplace or 
p-Laplace operators. Our aim is to investigate the cases, in terms of the equation parameters, 
for which the stated in the beginning property holds. 

We call the equation the model equation. It is is considered for x G (0,i?), < R< 00, 
with the conditions 

2/'(0) = 0, y{R) = 0, (2) 



and one looks for positive solutions. 
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We are interested in nonlinearities of power or exponential type. According to [7] we distin- 
guish the following two cases: 

- The Sobolev case: f{y) = Xy^, A = const, and 

a - /3 - 1 > 0; (3) 

- The Pohozaev-Trudinger case: f{y) = Ae^, A = const, and 

a = P + l. (4) 

If the equations involved the Laplace operator, both cases are related with the exceptional 
cases of the Sobolev theorem. In this regard we recall that the existence of solution to the 
equation Au + = depends essentially on the number p characterizing the nonlinearity. 
Namely, it admits positive solutions if and only if p > 2* = (n + 2)/(n — 2) — the critical 
Sobolev exponent [0]. 

For the Sobolev case the critical exponent associated to (1) was found in 0- Its value is 
given by 

_ (7 + l)(/3 + 2) 
^ - a-f3-l • 

Some explicit formulas for the solutions corresponding to the above cases were found in [Jj, 
among other results. The relations (jSJ and (jlj) were essentially used in that work. It is natural 
to raise the following questions: 

What is the relationship between the conditions © and Q and the nature of the equation (^? 
In what context do they appear? What is the role of the critical exponent? After all, why was 
it possible to find the exact solutions? 

We asked and answered these questions in |3 [3] . In those works the Lie point symmetry 
group of the model equation is calculated. Then it is verified which of the found symmetries are 
Noether symmetries. In this way, the property mentioned in the beginning regarding Noether 
symmetries and critical exponents, is established on the base of the S. Lie Theory. Since in 
the most cases the proof of embedding theorems between function spaces is reduced to radially 
symmetric functions, an ordinary differential equation is obtained. The symmetry approach 
permits as well the critical exponent be defined using directly that ordinary differential equation, 
without involving functional analysis. Namely: it is the only exponent for which a Lie point 
symmetry of Q with f{y) = Xy^ is a variational symmetry. 

Further, in |2] and [3] respectively, similar problems are considered for systems of ordinary 
differential equations and scalar partial differential equations. In this paper we shall present the 
main results of [H [3 IHl El • Our exposition in some parts follows closely the text of the original 
articles. It corresponds to the talk we gave on June 22, 2005, during the 6th International 
Conference 'Symmetry in Nonlinear Mathematical Physics', June 20-26 2005, Kyiv, Ukraine. 



2 The model equation 

The model equation contains the following particular cases: 

1. The equation is the radial form of partial differential equations containing 

- the Laplace operator in ifa = 7 = A^ — 1, /3 = 0; 

- the p-Laplace operator in ifa = 7 = A^ — 1, (3 = p — 2; 

- the /c-Hessian operator in iia = N — k, 'j = N — 1, P = k — 1. 



Noether Symmetries and Critical Exponents 



3 



2. The equation with a = ^ = N — l,(3 = and f{y) = Xe^ arises in the Liouville-Gelfand 
problem. 

3. The generalized Lane-Emden equation of the first kind: 

2/" + -y' + /3a:^~V = 0, (6) 

X 

can be obtained from ^ by setting P = 0, a = a, f{y) = (3y^, and 7 = a + — 1. The equation 
© with a = 2 is the Emden-Fowler equation, and if, in addition, = 1, then © is the equation 
proposed by Lane and studied in detail by Emden. 

This equation and its particular cases appear in astrophysics, mechanics, general relativity, 
various theories of gravitation, atomic physics and quantum mechanics. 

4. The generalized Lane-Emden equation of the second kind: 

y" + -y' + ^x'^'^e'^y = 0. 

X 

5. The Boltzmann equation 

m 

-{xy'Y = X ^ Xihe'^^y, 
1=1 

which is used in a huge number of applications, in particular, in the biophysics study of DNA 
modelled as a spiral on a cylinder . 



3 The Sobolev case for the model equation 

In this section we consider the model equation with f{y) = Xy^, A > and /3 > — 1, which 
can be written in the following form: 

^ /3+ 1 /5 + 1 ' ^ ^ ^ 

The symmetry group of ^ is calculated in [5], where various cases, depending on some 
relations between the parameters, are considered. Although that may be interesting from the 
point of view of the group analysis, we shall not present the complete results here directing the 
interested reader to [3] for further details. Rather, in order to establish the property stated in 
the introduction regarding Noether symmetries and critical exponents, we shall concentrate on 
the case when 

a - /3 - 1 > (8) 

and 

f3 + l<p<q*-l (9) 

where the critical exponent q* is given by ©. In this way, for the sake of clarity and brevity, we 
shall avoid some considerations which do not contribute to clarifying of that property. Moreover, 
these conditions are natural, since they are valid in the most of the interesting particular cases 
of the model equation. The following theorem is a direct corollary of the results obtained in 0. 

Theorem 1. Let (jH)) and © hold. Then the Lie point symmetry group of ^ is generated by 
the vector field 

p+i-p d d 

7 — a + p + 2 ox ay 
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Proof. By Q and© it follows that 7 — a + /3 + 2>0. (The latter condition is necessary for 
the existence of a positive solution of our problem 0.) It is also easy to see that Q and © 
imply a + 7 + [3^ 7^ 0. Then the theorem follows from the results in [3]. ■ 

Further we observe that the equation (jT)) is the Euler-Lagrange equation of the functional 
J[y\ = I L{x,y,y') dx, 



Jo 

where the function of Lagrange L = L{x, y, y') is given by 

L{x,y,y') = -^x^\y'f+^ — x^yP+\ 

Our next purpose is to check which of the symmetries leave the Euler functional J[y] invariant, 
that is, which of the already found in |Sj Lie point symmetries are actually variational. Again 
we consider only the case in the Theorem 1 in which all symmetries form a one-parameter group 
generated by Xi. 

Theorem 2 ( ) . A Lie point symmetry of the equation ^ is a variational symmetry if and 
only if p + 1 is equal to the critical exponent. 

Proof. We denote 

(3 + l-p 



k 



7-a + /3 + 2 

The Lie point transformation corresponding to (jlOj) is given by 

x* = fl'^X, 

y* = fJ'V- 

Clearly this is a scaling transformation. Then 

dy* 



dx*^ ' 



/3+2 ^ ^R* 



dx* — '— I (x*y(y*ix*)y^^dx* 



/3 + 2 Jo P+1 Jo 

Hence J[y*] = J[y] if and only if both exponents of /U vanish: 

/c(a-/3 -!)+/? + 2 = 0, 
/c(7 + l)+p + l = 0. 

The latter two equalities hold if and only ii p = q* — 1. ■ 

Henceforth, to the end of this section, we shall suppose p = q* — 1. Then by the Noether 
theorem we obtain that 

^ + ^ x^+'\y'f+^ - +1 + -^x'^^'y'^* = 0, (11) 



a-P-1 ' 7 + 1 

whenever y is a solution of the equation ((7| • 

Taking a look at [7], one can observe that the above first integral, up to some constant 
multiplier, is the i/j used in [71 p. 154]. Further we express 

_Xa;7-ay9*-i from ^ as a function 
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of X and y' and then substitute it into the equation ^ . We thus obtain that y' must satisfy the 
following equation 

" (7-0 + 1)1 , , 7 + 1 1 ,2 

y = — s— 1 y a — T~y ' 

p + l X a — p — 1 y 

whose solution, represented below, can be easily found by setting y' = vy. In this way we have 
come to the following theorem of [7j: 

Theorem 3 ([7j). The problem 

-{x'^\y'fy')' = Xx^y'i'~^ in (0,oo), 

y(0) = yo > 0, y'{0) = 0, (12) 
y > in [0, oo), 

has the unique solution 

y*{x) = {yo"" + kox')''^^" , 

where 

7-a + /? + 2 j-a + P + 2 ^ , P + l ( Xyo" V^^^^^^ 

a = , s = and kn = 

a-p-l ' /? + ! " a-/3-lV7 + ly 

The reason for the success in solving ()12[) is fact that the exponent in the right-hand side of 
the equation is critical, and hence any Lie point symmetry is variational, which reduces the order 
of integration procedure by two. In particular, this explains why the Lane-Emden equation 

9 

" I ' I 5 n 

y + -y + y =0 

X 

can be explicitly solved. Indeed, if the spatial dimension n = 3, then the corresponding equation 
parameters are 2 = n — 1 and 5 = (n -|- 2)/(n — 2) — the Sobolev exponent. 



4 The Pohozaev— Trudinger case for the model equation 

Denote by B the sphere in with centre at the origin and radius R. Then the Liouville-Gelfand 
problem consists of finding positive solutions of the following equation 

-Au = Ae" in B (13) 

with the Dirichlet boundary condition: 

u = on dB. (14) 

The equation above has been proposed and studied by Liouville. He found an exact solution in 
dimension = 1, and also (for = 2) a solution in terms of an arbitrary harmonic function. 
Bratu [2] found two explicit solutions of the problem if < A < 2/i?^ and N = 2. Gelfand [HI 
considered, among other things, the problem of thermal self-ignition of a chemically active 
mixture of gases in plane, cylindrical and spherical vessels. For iV = 3 he investigated the values 
of the parameter A for which the problem has a solution and studied the multiplicity of such 
solutions. 

It can be proved that if A < or if A is greater than a certain positive constant A* there is no 
solution of the boundary value problem (^SJi dJ)- Then by the celebrated results of Gidas, Ni 
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and Nirenberg JU] it follows that if a solution exists, then it must be radially symmetric. Thus 
the Liouville-Gelfand problem is reduced to the study of the following problem for an ordinary 
differential equation: 

-y" -E^y' = \ey in (0,i?), 

X 

y'(0) = y{R) = 0, 

y>0 in [0,i?). 
In [7j Clement, de Figueiredo and Mitidieri proposed the following generalization : 

-{x''\y'fy')' = \x^ey in (0,i?), 

y'(0) = y{R) = 0, (15) 

y>0 in [0,/?), 

assuming that 

a-/?- 1 = 0, /3>-l, 7>-l, 

that is, the situation corresponds to the Pohozaev-Trudinger case. 

Clement, de Figueiredo and Mitidieri found in [7j the above mentioned constant A*, proved 
that there is only one solution of (|15|1 if A = A*, and that there exist exactly two solutions if 
< A < A*. Moreover, using the method of first integrals, they found the explicit formulas for 
these solutions. Again the reason for the success of the integration is the fact that the condition 
a = [3 + 1 used in j7j holds if and only if all Lie point symmetries are variational symmetries, 
and hence the order of the integration procedure reduces by two. This is proved in where 
it is shown that the exact solutions can be accounted for by symmetry techniques. To see this, 
one needs first to find the Lie point symmetries of the model equation with f{y) = Ae", A > 0. 
We present the result in the following reduced form 

Theorem 4 {[^). Suppose that P > -I, a - P - 1 > and - a + P + 2 > 0. If P ^ or if 

P = 0, a ^ 1, the Lie point symmetry group of the equation 

-(x"|2/'|V)' = Ax^e^ 
is generated by 

1 d d 

X2 = X— + — , 

m ox ay 

where m = j — a + P + 2. If P = and a = 1 the symmetry group is the two-parameter Lie 
group determined by X2 and 

1/2 , \ d ^ d 

■ X — xmx — — h mx- 



7 + l\7+l J dx dy 

In the quasilinear case, it is easy to see by Theorem 4 that the following result holds. 
Theorem 5 (^y). Let P ^ 0. Then the Lie point symmetry group of 
-{x^Wfy'Y = Xx^ey 



is a variational symmetry group if and only if a = P + 1. 
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Further, observe that by Theorem 4, the Lie point symmetry group for (3 = and a = 1 
is two-dimensional. Hence, one can deduce that in this case there exists a one-dimensional 
subgroup consisting of variational symmetries. 

Then using a = /? + 1 and \y'\ = —y', we obtain the following first integral 

ipix,y,y') := ax"+^|yT+^ - (« + 1)(7 + 1)2:"|?/T + Ha + l)xT+^e^ = 0. (16) 

Again, taking a look at iTj, one can observe that the function ip above is the same that appears 
in [7j. Further we express y from ()16() as a function of x and y' and then substitute it into the 
equation in ()15|). We obtain that y' must satisfy the following Bernoulli equation 

„ (7-0 + 1)1 ,_l ,2 

y = y H —ry 

a X a + 1 
which can be easily solved. Thus 



y{x) = -{a + 1) In 



(a + l)(7 + l) 



+ C2. 



By the boundary condition y{R) = (see (2)) we obtain that 
ci = —j- and C2 = (a + Ij m ■ 



Aii(A^.)V" ^ ' (AM^)l/"' 

where fii, i = 1, 2, is a root of 

aAV<^ii(7+i)/o 



H{^x) := m^("+^)/° - + 1 = 0, m 



(a + l)(7 + l)("+i)/°' 
In this way we have come to another result of [7j: 

Theorem 6 ([7j). The solution(s) of (fT5|) can 6e represented in the following way: 

1) i/ < A < A* := (QT^iyj;7+i ^^en i/iere are two solutions 

2) if X = A* i/ien i/iere exists the unique solution 

Vo(a^) = —[a + 1) m \- — —r-, , , , x^ ^ " , 

^ \/^o (a + l)(°+i)/°i2(7+i)/° j 

where /Uq i^ the unique solution of H{fi) = 0; 

3) if X < or X > \* there is no solution. 

We observe that we have obtained the solutions found in [7] just by applying symmetry and 
variational methods, and without use of another first integral 

dx V 

See [7]. In fact, this ip corresponds to a dynamical symmetry. 
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5 The variational symmetries of the Lane— Emden systems 

In Sections 3 and 4, following |^ 1^, we have described the properties of the large class of 
quasilinear ordinary differential equations with power or exponential nonlinearities from the 
point of view of the Lie symmetry theory. We have shown that the Lie point symmetries of 
equations which contain critical exponents are actually variational symmetries. In this section 
we shall present the results of [HI which show that this property is also valid for systems. 
The system of two semilinear partial differential equations 

-Au = v'^, 

-Av = uP, (17) 

where the independent variable x € M", n > 3 — integer, p and q — positive real numbers, is 
called the Lane-Emden system. It can be considered as a natural extension of the celebrated 
Lane-Emden equation 

A9 + 9P = (18) 

in W. 

We are interested in radial ground state solutions of (fTTj) , whose study is reduced to the study 
of the following system of two ordinary differential equations: 

Tl — 1 

il H —ii + v'? = 0, 

V + ^^-j^v + = 0, (19) 



where t = \x\, u = etc., and n(0) = 'L'(O) = lim u{t) = lim v{t) = 0. 

t— +00 t—*oo 

It is well known that the Lane-Emden equation ()18() admits positive solutions if and only if 
p is greater than or equal to the Sobolev exponent (n + 2)/(?i — 2) (see, for instance, [0]). In [20] 
J. Serrin and H. Zou prove an analogous result for the system H17() . Namely: if {p,q) is on or 
above the curve in the first quadrant in the {p, q) plane given by 

+^ = n-2, (20) 



p+l q+l 

then there exist infinitely many (componentwise) positive radial solutions (u, v) of (|17|) tending 
to (0,0) when |x| — > oo. That is, the solutions are ground states. This theorem combined with 
the earlier nonexistence results of E. Mitidieri and J. Serrin, H. Zou ^1 gives a 

complete picture of existence and nonexistence of positive radial solutions of (|17jl in which the 
dividing curve (|20() plays an important role (see |201 Corollary 1.2, p. 370]). For this reason H20|) 
is called the critical hyperbola. 

In regard to existence of ground states for the Lane-Emden system H17|) . the results of 
P.-L. Lions J2] and J. Hulshof, R. van der Vorst imply that the problem with p and q 
satisfying ()20() has a unique, up to scalings and translations, ground state, which is positive, 
radially symmetric and decreasing in the radial component t = \x\. We shall come back to this 
point later. 

The basic assumptions on the parameters are: 

n>3, pq>l, p^q, p ^ 1, q ^ 1. (21) 



Then the symmetry group of the Lane-Emden system (|19|) is calculated in [H], and we state it 
in the form of the following 
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Theorem 7 (0). Suppose that the conditions ()21() hold. Then the Lie point symmetry group 
of the Lane-Emden system l|19jl is the one-parameter group generated by the vector field 

^ d 2(1 + q) d 2{l+p) d 
at I — pq ou 1 — pq ov 

Further, following the same argument as in the proof of Theorem 2, we obtain the next 

Theorem 8 (16J). Let the parameters n, p and q satisfy ()21() . Then any Lie point symmetry 
of the Lane-Emden system (|19j) is a variational symmetry if and only if 

2 pq-l ' 

that is, the point {p, q) is on the critical hyperbola l!\2(J\i . 

The statement of Theorem 8 is a conjecture of Enzo Mitidieri |16j . 

We observe that if {p, q) belongs to the critical hyperbola, the Noether Theorem immediately 
gives the first integral used in and called invariant paraboloid. 

We conclude this section noting that the exact solution of the problem ()19|) is not available. 
As mentioned above, it is known to exist, and only the asymptotic behavior of such ground state 
is known. Challenged by Enzo Mitidieri and provoked by J. Hulshof, R. van der Vorst's paper 
we proposed in [2] the following 

Problem. Find the solution of the Lane-Emden system explicitly. 

6 The Lie point symmetries of semilinear polyharmonic 
equations 

In the next section we shall consider the semilinear polyharmonic equation 

(-l)'=A^ = /(n), (22) 



where A is the Laplace operator in M", n > 2 and k > 1. The Lie point symmetries of ((22} 
were studied by Svirshchevskii, [^HIISI who obtained the complete group classification. We shall 
cite partially his result in a form suitable for our purpose to investigate the cases of power and 
exponential nonlinearities. 

Theorem 9 ( \22\ I23| ). The widest Lie point symmetry group admitted by (|22|) with general 
f{u) is determined by translations and rotations. 

For some special choices of the right-hand side f{u) it can be expanded by additional operators 
as follows: 

U fiu) = u^ , p ^ 0, p I, we have the generator of dilations 

d 2k d 

^ = 

OXi 1 — p ou 

and, for p = (n + 2k) /{n — 2k), there are n additional generators 

d d 

Yi = {2xiXj - \x\'^6ij) — + {2k - n) Xi u—. 

Further, if f{u) = e" then the operator 

d d 
W = Xi — -2k — 
OXi ou 
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generates a sub-group of the Lie point symmetry group of <\2'^. If n = 2k, there are n additional 
generators 

d d 
Vi = {2xiXj - \x\'^6ij) — - Akxi —. 

Above is the Kronecker symbol, = ( ^^2, ^s] ^'^d summation over a repeated index 
is assumed. 



7 The Noether symmetries of semilinear polyharmonic 
equations with critical nonhnearities 

Up to this point, we have estabhshed the property relating the critical exponents and the Noether 
symmetries for some ordinary differential equations and systems. It is natural to conjecture that 
the Lie point symmetries of more general differential equations involving critical exponents are 
Noether symmetries. In |2] we show that this is valid for a class of partially differential equations. 
Namely, for the semilinear polyharmonic equation (|22() . 
Again, we shall distinguish two cases: we have 

/(n) = tiP, n>2k>2, p^O, p^l, (23) 

in the Sobolev case, while in the Pohozaev-Trudinger case 

f{u) = e", n = 2k, k>2. 

Having at our disposal the already cited group classification (see the preceding section), 
we investigate which of the Lie point symmetries in these cases are variational or divergence 
symmetries. 

The first results in can be formulated as follows: 

Theorem 10 (P])- Suppose that the conditions l\2'6^ hold. Then the Lie point symmetry Z (see 
Theorem 9) of the polyharmonic equation 

{-!)'' A'^u = 

is a variational symmetry if and only if 

n + 2k 
n — 2k 

that is, p is equal to the well known critical exponent. 

Theorem 11 {^). Let n > 3 and k >1. Then the Lie point symmetry W (see Theorem 9) of 
the polyharmonic equation 

(-l)'=A'=n = e". 

is a variational symmetry if and only if 



n = 2k. 
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The importance of the relation n = 2k for the polyharmonic equation was pointed out to 
us by Enzo Mitidieri in January 1995 jl4j . 

The proofs of Theorems 10 and 11 are analogous to the proof of Theorem 2. In fact, this 
is a dimensional analysis argument since the latter is equivalent to the invariance under scaling 
transformations (see, for example, [Tl [T7]^. 

Since the translations and the rotations are variational symmetries of the considered polyhar- 
monic equation, it follows from Theorems 10 and 11 that in the critical cases the corresponding 
Lie point symmetry groups contain a subgroup consisting of variational symmetries. However, 
in the critical cases one can get more. 

The main results in P] are the following 

Theorem 12 ([3|). Suppose that the conditions ()23() hold. Then any Lie point symmetry of the 
polyharmonic equation 



is a divergence symmetry. 

Theorem 13 (jSl). Let n > A, k > 2 and n = 

polyharmonic equation 

(-l)'=A'=n = e" 



2k. Then any Lie point symmetry of the 



is a divergence symmetry. 

The proofs of Theorems 12 and 13 are reduced to the proof of the facts that Yi and Vi (see 
Theorem 9) are divergence symmetries, since the other infinitesimal generators of the corre- 
sponding Lie point symmetry groups determine variational symmetries by Theorems 10 and 11 
as already observed. For this purpose we find explicitly the vector-valued 'potential' function B 
determining Yi and Vi as divergence symmetries [3]. 

The importance of the variational and the divergence symmetries is due to the fact that 
they determine conservation laws via the Noether Theorem ^ El • Thus the next step in this 
research is to establish the conservation laws corresponding to the already studied variational 
and divergence symmetries. This is possible since we already have at our disposal the explicit 
formula for B. This is also done in [21 • 

The Theorems 10-13 confirm the validity of the general property relating Noether symmetries 
and the critical parameters, stated in the beginning of the introduction. Another examples which 
illustrate this property are given in in which Svirshchevskii proved that the symmetries of 
the p-Laplace equation 

-divdV-ulP-^Vn) = \u\''~'^u 
and the equation 

-A(|An['"An) = jul^-^-u 
are variational if and only if the parameters assume critical values. 
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